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Abstract. The standard model of cosmology is based on homogeneous-isotropic 
solutions of Einstein's equations. These solutions are known to be gravitationally 
unstable to local inhomogeneous perturbations, commonly described as evolving on 
a background given by the same solutions. In this picture, the FLRW backgrounds 
are taken to describe the average over inhomogeneous perturbations for all times. We 
study in the present article the (in)stability of FLRW dust backgrounds within a class 
of averaged inhomogeneous cosmologies. We examine the phase portraits of the latter, 
discuss their fixed points and orbital structure and provide detailed illustrations. We 
show that FLRW cosmologies are unstable in some relevant cases: averaged models are 
driven away from them through structure formation and accelerated expansion. We 
find support for the proposal that the dark components of the FLRW framework may 
be associated to these instability sectors. Our conclusion is that FLRW cosmologies 
have to be considered critically for their role to serve as reliable models for the physical 
background. 



PACS numbers: 04.20.-q, 04.25.-D, 04.40.-b, 95.35. +d, 95.36,+x, 98.80.Es, 98.80.Cq 
1. The problem of a physical background in cosmology 

The standard description of the global evolution of the universe relies on the class 
of homogeneous-isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions 
of Einstein's equations. By construction all inhomogeneities in matter and geometry, 
responsible for the formation of structures, average out on this given background. The 
conjecture is held that the FLRW background is the actual physical background. 

This widely adopted framework needs to be challenged for several reasons. First, 
apart from the question of how to technically implement an average (a non-trivial 
subject for tensors), the result will depend in any case on the spatial domain over 
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which inhomogeneities are considered: a scale dependence of the averaged variables 
must generically result, casting strong doubts on the scale-independent averaged values 
issued from the FLRW solutions (see [20] for a summary of these thoughts). Secondly, 
it is actually naive to expect that a strictly homogeneous model provides the average, 
given the non-linearity of Einstein's equations and its implications (e.g. the non- 
commutativity of the time-evolution and averaging). This background issue may also 
be identified as the crucial question in the discussion of whether or not the need for 
dark components can be replaced by employing inhomogeneous models |22j . 

Looking at current cosmological structure formation models based on perturbation 
theories or iV-body simulations, we find that the extra (backreaction) terms due to the 
averaging procedure vanish on the background of a homogeneous-isotropic solution. A 
closer look reveals that this property results by construction rather than by derivation, 
and it is a consequence of employing Newtonian or quasi-Newtonian schemes with 
periodic boundary conditions and a flat background geometry (see [9] for the proof). 
This construction is not expected to work in the framework of general relativity because 
of: (i) the relevance of the spatial intrinsic curvature (the second derivatives of the 
metric may be significant even if the metric perturbations are negligible [ID]), together 
with (ii) the fact that inhomogeneities are coupled to the spatial curvature evolution [3], 
and finally (iii) the absence of a conservation law for the averaged intrinsic curvature [8]. 

It is of course conceivable that the homogeneous solutions provide in some 
spatial and temporal regimes a good approximation for the evolution of the averaged 
distribution. A systematic approach to address this assumption is to analyze the 
stability of the FLRW solutions in the space of averaged cosmologies. Although FLRW 
solutions are known to be unstable to local perturbations, we may be able to find that 
they enjoy stability if subjected to perturbations on average and on some large scale. 
In the present paper we follow this approach, and we find that the FLRW solutions 
can also be globally unstable within a class of averaged models. Structure formation 
and accelerated expansion imply that the physical background is driven away from the 
FLRW background in relevant cases. Our investigation is independent of the choice of 
scale, and it therefore contains all the dynamical situations a domain can undergo. 

For the class of inhomogeneous cosmologies we consider in this work, we shall 
employ a framework that averages the scalar parts of Einstein's equations, in which 
case the averaged model is unambiguously defined (reviewed in [6j[7]) and, due to its 
covariant implementation, is under control pQEUEHEE]- ^ e are §°i n § to explore the 
phase space of averaged cosmologies by working within a class of exact scaling solutions, 
generalizing previous results j . 

The study of the phase space of dynamical systems has a long tradition. The 

t I n [11] - the scaling law was applied to the total averaged curvature, and wc choose here to apply it 
to the averaged curvature deviation from a constant (Friedmann) curvature. Taking a null constant 
curvature in our work leads back to the results of [llj . This difference of approach allows us to study 
in this paper the stability of all FLRW backgrounds, not fully addressed in [TTJ. Note already that the 
phase portraits drawn in [TT] correspond to the one-dimensional locus Stf? = in the plots of figure [1] 
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most comprehensive study of FLRW cosmologies that include radiation, dust and a 
cosmological constant has been provided by Ehlers and Rindler [19]. Other studies, 
including Bianchi models, were treated by Wainwright and Ellis [28] who proposed the 
so-called dynamical system approach, which we shall use in this paper. The same 
technique can be followed to investigate the dynamics of minimally coupled scalar fields 
in the FLRW framework [18]. Note that the backreaction fluid we study in this paper 
can also be interpreted as an effective scalar field [11] that allows comparison of this 
work with our analysis. More recently, dynamical system analysis was used for modified 
theories of gravity (fourth order [13111411161123] . scalar tensor [15] and Hofava-Lifshitz |17j ) 
leading to a greater understanding of the cosmology of these theories. Also, results of a 
recent investigation of LTB models [27] are in accord with our finding^. 

In the next section, we introduce the basic framework of averaged inhomogeneous 
cosmologies, choosing a vanishing cosmological constant for simplicity. We consider 
in section [3] the backreaction terms to obey scaling evolution laws, and present and 
analyze the corresponding autonomous dynamical system. In section HI we summarize 
and discuss the results, and we conclude in section [5j 



2. Effective description of inhomogeneous cosmologies 

Let us consider an inhomogeneous, irrotational and pressureless fluid (dust). The 
spatial average over a compact, restmass preserving domain T> of the Raychaudhuri 
equation, the Hamilton constraint and the dust continuity equation reads in its rest 
frame (see (3jlU|6] for details) 

a v AtiG Qv 

— + -5- (Q)v = -5~> W 
ax> o o 

a v \ 2 8irG k^ Wy + Qv , , 



Cl£> J 6 Ojy 

{Q)v + Z — (q)v = °> (3) 
ax> 

where the angular brackets stand for the Riemannian spatial volume average, and the 
overdot is the partial time derivative, here identical to the covariant derivative, av is 
the effective scale factor, (g) v is the averaged energy density of the dust, Qv is the 
kinematical backreaction term, and Wv is the averaged curvature deviation from a 
constant curvature, defined as 

av ■= (Vv/V Vi ) 1/3 , (q) v = (g) Vi a v 3 , (4) 

Qv ■= I ((9 - (6} v ) 2 ) v - 2 (a 2 ) v , W v := (R) v - Qk Vi a v 2 , (5) 

with Vv being the volume of the domain, TZ its three-Ricci scalar curvature and 9 and 
a the expansion and shear rates of the dust. (Here and in the following, the subscript 

§ In [57], Sussman used quasi-local variables instead of averaged ones. The relation he gave between 
them allows us to establish common findings. 
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% stands for the initial value.) Equation (j2J) is an integral of equation (CQ) under the 
conservation law Q, if and only if Qv and Wv obey the integrability condition 



The set of equations (p3l3]) formally resembles, from a kinematical point of view, a 
Friedmann cosmology sourced by two fluids, but there are some fundamental differences 
between the two frameworks. First, the system is valid for any metric, whereas the 
FLRW equations only hold for maximally symmetric metrics. Second, the second fluid 
(the backreaction fluid) emerges from the averaging process; it is not added in our model 
contrary to a two-fluid FLRW cosmology. This fluid represents the difference between 
the dynamics of a Friedmann dust background and of a physical dust background 
obtained through the averaging proceduretJJ]. Third, the dynamics of the averaged 
geometrical variables of our model generally differs from that of the FLRW model: 
(1Z) v might evolve differently from the Friedmann curvature kvfi^ restricted to the 
same domain. Fourth and last, the averaging procedure is performed over a chosen 
domain, and it implies a scale dependence of all the variables. In particular, one expects 
the backreaction fluid to exhibit rather distinct behaviors according to the scale |29j . 

For void domains (g — 0) the backreaction terms, encoded through geometrical 
invariants, correspond to inhomogeneities of the spatial geometry of the vacuum. For 
domains over which Qv = 0, inhomogeneities could still be locally present even though 
the kinematical backreaction identically vanishes on T> due to an exact compensation 
between the expansion fluctuations and the sheajjjl In the same spirit Wv = only 
states that the averaged curvature deviation vanishes on D, but nothing can be inferred 
about the local curvature deviations within the domain. Homogeneity is only achieved, 
if the backreaction variables are required to vanish on all scales. 

Since we aim at studying the instability of Friedmann dust backgrounds, we shall 
consider for our analysis the quantity X v ■— Qv + VVp, which corresponds to the whole 
averaged departure from the Friedmann framework. From the constraint (El), a null X-p 
implies vanishing Qv and Wp, in which case we recover a FLRW background on average. 
Also, the dark components can be thought of as the manifestations of the deviation of 
the physical background from a FLRW one: a positive Qv contributes to accelerate the 
expansion of the domain and plays against gravity, whereas a negative Qv participates 
in the deceleration of the domain's expansion and adds to gravity (see equation (TjQ)). In 
this sense, the instability sectors, defined through Xv, represent the dark components 
of the concordance model. 

|| This is an essential remark when one considers perturbations, the choice of the background being 
crucial. 

For example, this happens for zero curvature, spherically symmetric Lcmaitrc-Tolman-Bondi 
solutions, see [7] in this volume. 




{Qv a v ) ' + a v 2 



(WvalY =0. 



(6) 
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3. Dynamical system for the scaling solutions 

3.1. The scaling solutions 

We need one additional relation to solve the system ([T]l3]), and this closure relation 
expresses the freedom of the choice for the local structure of the inhomogeneities, 
encoded on average by the dynamics of the backreaction fluid3- A natural choice is 
given by establishing a formal analogy between the backreaction terms and a fluid, 
whose energy density and pressure depend on the volume of the domain [TT] : 

QvocVS, Wp oc => Q v = Q Vi al , W v = W Vi a p v , (7) 

where the scaling parameters n and p are real constants. Expressions (ED) together with 
([6]) result in the following constraint the dynamical system has to satisfy: 

(n + 6) Q Vl a n v + (p + 2) W Vi a p v = . (8) 

In what follows we shall concentrate on the class of scaling solutions n = £@, for which 
relation (JBj) implies 

{n + 2)W v = -(n + 6)Q v (n + 2)X V = -AQ V . (9) 



3.2. The autonomous system 

Let us introduce the Hubble functional Hx> '■= av/av an d the volume deceleration 
functional qx> ■= — 1 — Hu/H v . We define domain-dependent dimensionless cosmological 
characteristics as 



oil j) a-p ii-p vii'D 

which by construction add up to 1 according to the Hamilton constraint (J2J). Upon 
using the expressions (MTU]) we replace the system (fI U5fB]) by 

(11) 
(12) 
(13) 
(14) 

where the prime denotes the derivative with respect to the evolution parameter 
Nd '■= la.au- The two last equations are easily derived noting that (1/H V )' = 2(qu + 
1)/H V . The autonomous system f lT3|T4|) determines the orbit of a cosmological state, 
defined by the quartet (fi^, &x, Of , n), in the corresponding phase space, and its fixed 
points are provided in table [U The phase space (f2^, fi^) presents three invariant lines 
under the phase flow: ti^ = 0, fi£ = and fi£ = 1 - fi£ - fig = 0. 

+ The need for a closure relation is also present in the standard FLRW framework, and it corresponds 
to the choice of an equation of state for the fluid sources. 

* Note that the unique solution for n ^ p (n = — 6,p = —2) is physically equivalent to the case 
n = p = —6: Qx> follows the same evolution law in both situations, and so does the total averaged 
curvature (72.) j, which evolves proportionally to a^, 2 . 
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Table 1. Fixed points of the dynamical system (|13|14j) . The sign of their eigenvalues 
determines their stability pp. 





Coordinates 


Scale factor 






Stability 

{eigenvalues} 












n < — 2, attractor 


point A 


(0,0) 


a v (t) = H Vi (t-U) - 


r 1 


2 


n > —2, saddle 

{-1,71 + 2} 

n < —3, saddle 


point $ 


(1 0) 


a-n(t) = (iHn (t - U) - 


r 1 

_2 


.3 


Tt > — 3, repeller 
{l,n + 3} 

?i < —3, repeller 


point C 


(0,1) 


a v {t) = {-^H Vl (t-U) + l)~ 
a v (t) = cxp(H Vt (t-U)) 


11 

(" 


(n ± 0) 
= 0) 


—3 < ?i < —2, saddle 
71 > —2, attractor 
{-77 - 3, -?i - 2} 


line Ci 




a v (t) = (t - U) - 


HI) 


2/3 


(77 = -3) 


line £2 


^ = 


a v (t) = H Vi (t-U) 


+ 1 




(n = -2) 



There are five situations for which the stability of the fixed points does not change: 
n < —3, n = —3, —3 < n < —2, n = —2 and n > —2. The fixed point A shares 
the same properties as a Milne domain (i.e. a void domain with FLRW kinematics and 
FLRW curvature on average); the fixed point B corresponds to an Einstein-de Sitter 
dust domain (i.e. FLRW kinematics and zero FLRW curvature), and the fixed point C 
corresponds to domains filled only with the backreaction fluid. The stability of A and B 
can change according to the value of the scaling parameter, but they always depict the 
same cosmological models. In return, C always describes different cosmological models 
which might also have a different stability. For the bifurcation value n = —3, there is a 
transfer of stability between B and C which both belong to the line of fixed points £1, 
while at the bifurcation value n = —2, there is a transfer of stability between A and B, 
both belonging to the line £2- 

Finally, note that the dynamical system (I13|I14I) is written according to the evolution 
parameter Nx> and not according to a time parameter. Consequently, the phase portraits 
drawn from this system do not suffice to deduce the time evolution of a cosmological 
state, since we also need information about the time evolution of Nx> or, equivalently, 
ax>- We shall therefore also indicate for all phase portraits the behavior of the scale 
factor. 
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3.3. Phase portraits 

The phase portraits of the five cases are shown in figured], and the numbers (background 
colors) of their regions correspond to the evolutions of ap given in figure EJ To deduce 
the time evolution of a state, one shall refer to these two sets of plots. For the sake of 
clarity, we provide hereafter a detailed explanation of the correspondence between them 
in the situation n < —3 (first plot of figured]). 

A state &x) belonging to the region I (red) evolves from the fixed point C when 
its scale factor is zero (plot I (red) of figure [2]), through its initial conditions, to infinity 
when its scale factor attains its maximum. This divergence is due to the singularity 
of the cosmological parameters for Hp = and not due to a physical boundless blow of 
the domain. Having reached its maximal extension, the domain shrinks and (f2^,f2^) 
comes back from infinity and tends to C, while ap decreases from its maximum to zero. 
Hence, the fixed point C is both a past repeller and a future attractor. 

In the region II. a (green), the states evolve from C, when ap is zero, toward A, 
when ap tends to infinity (plot II. a (green) of figure E}, without any accelerated period. 
In this region, C is a past repeller and A is a future attractor. 

A state in IV (blue) leaves the fixed point A, while the scale factor decreases from 
infinity to its minimum in the past, and reaches infinity. It then goes back to A when 
the scale factor grows to infinity (plot IV. a (blue) of figure [2]). The state undergoes an 
accelerated period when operating its bouncing back from infinity. Here, A is a past 
repeller and a future attractor. 

Finally, in the region V (brown), (fi^> fiy) runs from infinity at the minimum of ap 
to infinity at its maximum, and it oscillates (plot V (brown) of figure [2]). The minimum 
of ap is reached when the orbit runs below the dashed line since it corresponds to a 
positive acceleration. 

All the orbits that go to infinity only describe part of the evolution, i.e. up to a 
maximum or minimum of the scale factor. The rest of the evolution is described not 
by reversing the direction of time, but by reversing the monotonicity of ap. In other 
words, the state travels back on the orbit of the phase portraits. 

The other generic orbits expounded in figured] are the following (see also figure |2]): 
a state in IV.b (blue) follows the same type of evolution as a state in IV. a except that 
it always undergoes an acceleration; the evolution of a state in II. b (green) differs from 
that of II. a by the fact that it experiences an acceleration from a finite time; and a 
state in III (darkest green) evolves from B, when ap is zero, toward C, when ap tends 
to infinity. In this region, B is a past repeller and C is a future attractor. 

Finally, we present the evolution of the scale factor for domains lying on an invariant 
line. If a state belongs to the segment between two fixed points, it experiences a big- 
bang and a future expanding phase. If it belongs to a semifine, either the domain 
undergoes a big-bang and big-crunch evolution when the arrow on the invariant line 
leaves the fixed point, or the domain undergoes a past and future expanding phase when 
the arrow is orientated toward the fixed point. 
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Figure 1. Phase space (^„, &x) of the dynamical system (|13ll4j) . From left to right 
and top to bottom: n < —3 (n = —4 for the illustration), n = —3, —3 < n < —2 
(n = —2.5), n = —2 and n > —2 (n = —1). The thick straight lines (dark blue) arc 
the invariant lines fl^ = 0, Q^- = and fl^ = 1 — — ^ = 0. Every line parallel 
to the latter, the inclined line, corresponds to a constant with values increasing 
downward. The dashed line shows a vanishing qx>] it is an orbit only for n ~ — 3, and 
it coincides with — for n = — 2. A domain is accelerated below this line and 
decelerated above when n < —2, and this correspondence is inverted for n > — 2. The 
arrows show an increasing TV-p , the dots are the fixed points of table [TJ and the square 
in the last plot stands for the observational values (see section [473]) . 
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Figure 2. Evolution of the scale factor of a domain with respect to Hx> t t. The numbers 
(colors) refer to the different regions of the phase spaces in figure [1] All the states of 
a region have a scale factor which evolves in the same way. We have chosen a positive 
Hu i to deal with initially expanding domains and a null U. Region I (red) shows big- 
bang and big-crunch domains, regions II show big-bang and future expanding (with 
an asymptotically null acceleration) domains, region III. 1 (darkest green) describes 
big-bang and future expanding (accelerating) domains for —2 < n < 0, region III. 2 
(darkest green) describes big-bang and big-rip domains for n > 0, regions IV (blue) 
show past and future expanding domains, and region V (brown) shows oscillating 
domains. We distinguish two subcases for the regions II and IV: for the former, region 
II. a is never accelerated, while region II. b is accelerated from a finite time onward; for 
the latter, region IV. a is accelerated during a finite period of time, while region IV. b 
is always accelerated. 
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4. Summary and discussion 

First and foremost, we clarify the definitions of the different types of domains at stake. 
A domain T> is Friedmann, i.e. locally homogeneous and isotropic, if and only if for 
all subdomains £ C T> one has Qq = f2yy = or, equivalently, VL X = 0. We have 
focused our analysis on one single domain T> without any multiscale considerations; so 
it is fair to assert that the Friedmann realization is not achieved in our approach since 
no additional suppositions are made. A domain with a null Q x and with nothing more 
inferred about the Q x of the subdomains, is most probably non-Friedmann: as noted 
in section [2} the cancellation of the backreaction fluid can be the result of an exact 
compensation between the expansion variance and the shear, and inhomogeneities and 
anisotropies may still be locally present. Nevertheless, such domains are destined to 
follow on average a Friedmann evolution since their inhomogeneities do not contribute 
on average to the kinematics. We shall call these domains Friedmann-like domains. 
Therefore, we shall deal with Milne-like and Einstein-de Sitter-like domains, instead of 
locally Milne and Einstein-de Sitter domains. This complies with the actual physical 
situation: Friedmann models as dynamical models for the average evolution are only 
realized for strictly homogeneous distributions of matter and geometry, but there exist 
situations — called Friedmann-like — in which Friedmann models do provide the evolution 
on average, i.e. in which they describe the physical background. 

Another terminology we employ makes use of equation ([9]): for n < —2, Q x > 
corresponds to shear-dominated domains, and Q x < labels expansion variance- 
dominated domains. This correspondence is inverted for n > — 2. 

4-1. Summary of results 

Before delving into details, simply recall that the orbits strongly depend on the leading 
term in the cosmic trio oc a^, 3 if^ 2 , oc a^H^ 2 and Q x oc a^H^ 2 . When 
n < —2, the constant curvature kv t , useful to compare the effective background with 
the Friedmann-like background, determines the dynamics of the domain, while in the 
remaining situations the backreaction fluid is the relevant quantity. 

For n < —3 and for a negative constant-curvature fcx> s , shear-dominated and 
expansion variance-dominated domains are attracted by a Milne-like state. In the past, 
the former originates from a state with Q x ~ 1, whereas the latter comes from a Milne- 
like state. Expansion variance-dominated domains are therefore asymptotically, in the 
past and in the future, dominated by a negative constant curvature. The bounce of these 
domains occurs when the contribution of the backreaction fluid is no longer negligible. 
For a positive k-p v expansion variance-dominated domains are oscillating (this is the 
only situation where it occurs): when the volume of the domain is sufficiently small, the 
backreaction fluid — mimicking a dark energy behavior — becomes preponderant over the 
other components, thus avoiding the collapse. In contrast, shear-dominated domains — 
mimicking a dark matter behavior — evolve toward a domain filled only with Q x and 
collapse in a finite time. 
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At n — —3, all domains with a negative fcx> s approach a Milne-like state, and all 
domains with a positive kx> % eventually collapse. The attractor state of these latter 
depends on the initial conditions: though it always exhibits f2f = 0, all the values 
of fi^ and 0% are attainable provided the Hamilton constraint is satisfied. For the 
particular state Q^ 1 = 1 (dashed line of the second plot of figure EJ), the domain 
experiences a stationary expansion (d© = const). always being equal to 1 in this 
situation, the dust and backreaction fluid energies exactly compensate each other along 
the whole evolution. Furthermore, we here confirm that this stationary state is stable: 
any perturbation always leads back to a stationary state (as was conjectured in [5], see 
also [2}|25] for comparisons). 

For —3 < n < —2, shear-dominated and expansion variance-dominated domains 
have the same qualitative asymptotic properties. For a negative kx>i, they originate from 
an Einstein-de Sitter-like dust state and are attracted by a Milne-like state, whereas for 
a positive constant curvature they emanate from and converge to an Einstein-de Sitter- 
like state and eventually collapse. In these last scenarios, even if the backreaction fluid 
mimics a dark energy behavior over the domain, its intensity is not sufficient enough to 
counterbalance the collapse contrary to the case n < — 3. 

In all the previous situations (except for the region V (brown)), a positive constant- 
curvature domain is bound to collapse, and a negative constant-curvature domain always 
converges toward a Milne-like state. Therefore, a small perturbation around a null kv, 
seals a totally different fate according to its sign. 

For n = —2, domains with Q^' < 1 are attracted by a Milne-like state, and domains 
with > 1 converge to an Einstein-de Sitter-like state and collapse within a finite 
time. The constraint ([6]) here implies Qx> = and oc a^> 2 , which means that there is 
no kinematical backreaction, but also no curvature deviation from a constant curvature 
stricto sensu. All states follow the kinematics of a Friedmann-like state on average, and 
the whole drawn phase portrait is physically equivalent to the line = 0. 

Finally, for n > —2, shear-dominated domains emerge from and are attracted by 
an Einstein-de Sitter-like dust state, and they eventually collapse. Expansion variance- 
dominated domains also originate from an Einstein-de Sitter-like state, but they are 
attracted by a state filled only with the backreaction fluid, and they experience an 
accelerating expansion from a finite time onward. As an example for the latter: void 
domains that are nearly Friedmann on average (f2^ = 0, > 0) are unstable and 
driven away toward this expanding state; this instability is the origin of the possible 
onset of an inflationary scenario created from initial curvature inhomogeneities of the 
Einstein vacuum [12] . 
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4-2. Discussion: instability sectors 

A state leaves the class of FLRW backgrounds and enters into what we define as an 
instability secto^ as soon as the backreaction fluid contributes to the kinematics of the 
domain. The existence of a non-vanishing averaged contribution of inhomogeneities 
invokes deviations of a cosmic state from its Friedmann-like characteristics. Given the 
generic behavior of any domain, we now address the issue of the stability of the FLRW 
backgrounds: does a Friedmann-like state, lying on the line fi^c = 0> converge on average 
to the same or another Friedmann-like state when subjected to perturbations? 

A convergence toward any Friedmann-like state or, equivalently, toward the class of 
FLRW backgrounds, happens when the two following conditions are met: first, X-p — > 0, 
which entails that inhomogeneities are negligible on average; second, — > 0, which 
states that the energy contribution of the backreaction fluid becomes subdominant on 
average compared to the other components. Note that for the particular case n = —2, 
as explained above, the averaged contribution of inhomogeneities vanishes, and all the 
states are Friedmann-like. This situation does not admit any global inhomogeneity 
effect, and we shall therefore disregard it in our discussion. 

Based on the previous analyses, we conclude that FLRW backgrounds are stable 
for: (i) n < —2 and fcp. < 0, the Milne-like state is an attractor for any inhomogeneous 
deviations, and (ii) n < —3 and kj) i = 0, the Einstein-de Sitter-like dust state is 
eventually reached. In these situations the FLRW backgrounds are stable, in the 
situation (i) the background can only be attained asymptotically. 

In the other situations, the perturbations of any Friedmann-like state lead either 
to Q^- 4 1 or Xp — > oo or both, namely to a non Friedmann-like state. The 
FLRW backgrounds are therefore globally unstable for: (iii) n < —2 and kv i > 0, 
(iv) —3 < n < —2 and kp> i = 0, and (v) n > —2. In these situations they obviously 
constitute an incorrect approximation to the physical background. 

4-3. Discussion: interpretation as dark sectors 

The different cosmological observations interpreted in the context of FLRW cosmologies 
constrain the instability sectors to which the dark components could be associated. For 
instance, we may put a conservative upper bound on the present-day cosmic state, 
^m,toda y < 1> i- e - an y physical orbit should respect this bound, if we consider small- 
scale (galaxy clusters and voids) or large-scale (CMB and high-redshift supernovae) 
cosmological domains. For instance, this restriction forbids some of the orbits for these 
gravitational systems: they cannot belong to the region V (brown) of figure [U and so 
we cannot have oscillating domains. 

The dark matter component can be associated with domains that are dominated by 

jj From figure [1] it is obvious that no Friedmann-like state is stable (except the point A), i.e. in those 
unstable cases no perturbation slightly above the fip = line relaxes to its original point of the phase 
space diagram. Therefore, we address the more general issue of stability of the class of backgrounds 
with Friedmann-like behavior. 
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shear fluctuations, and which eventually collapse. The instability sectors of the regions 
I (red) in the first plot of figure [fl in the second and third plots (with VL^ > 0), and 
in the fifth plot (with fi^c < 0) may present an effective dark matter behavior for the 
backreaction fluid. Concerning the dark energy component, it is effectively mimicked in 
the instability sectors of the region IV. b (blue) in the second plot, II. b (green) in the 
third plot, and III (darkest green) in the last plot. In each of these cases, there is a global 
acceleration of the domain. Note, interestingly, that the sign of the third derivative of 
a© may also change, allowing for a variety of different dynamical evolution histories 
of the dark energy component. In the following, we pick two proto-type examples of 
instability sectors to which the dark components can effectively be associated. 

Consider a small domain filled mainly with dust (f2^ ! ~ 1), dominated by its 
shear fluctuation (Q 1 ^ > 0), with an initially slightly positive constant curvature 
(Q^ 1 < 0) and such that n < — 3. This situation can correspond to negative or 
positive averaged scalar curvature (TZ)^, and it describes a pancake-like configuration 
of dust, like a proto-sheet (or proto-filament). Such a domain evolves toward a 
state dominated by the backreaction fluid, passing close to the observational ratio 
^m/^x = ^baryon/^DM — 0.21, and eventually collapse iffl (see figures [T] and . This 
example shows that the backreaction terms can turn out to play the role of a dark 
matter component, thereby accelerating the formation of structures by adding gravity 
and accreting the dust. 

Consider a large domain, say our Hubble sphere, endowed with an initially slightly 
negative constant curvature (fi^ 1 > 0) and mainly composed of dust ~ 1). Suppose 
also that the expansion variance slowly varies and dominates the backreaction fluid (so 
—2 < n < and > 0). Such a domain is always characterized by a negative averaged 
scalar curvature {TZ) V and enters, from a finite time onward, a phase of accelerated 
expansion, may pass close to the observational value (fi^, — (0-27, 0.73) (the square 
in figure H]), and eventually ends up being entirely dominated by the backreaction fluid. 
Recall that the backreaction fluid can be formulated as a scalar field [11], and here 
its equation of state would be given by u = —l — n/3 > —1. These accelerated 
expansion scenarios include the case found in perturbation theory [23], and they would 
be interpreted as the signature of dark energy in the standard model. 

These two examples, though not depicting exhaustively all the possible scenarios, 
show how the instability sectors of FLRW backgrounds can be linked to the dark 
components. 

5. Concluding remarks 

We have studied in the present paper the global gravitational stability of FLRW 
backgrounds in the class of scaling averaged inhomogeneous models. Our investigation 
has led to a detailed classification of stability and instability sectors in the phase space 

ft Note that our effective model does not take into consideration either the pressure or the vorticity of 
the matter at stake, which obviously cannot be neglected at the end of the collapse process. 
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of the scaling averaged cosmologies, and it has shown that averaged models are driven 
away from FLRW backgrounds through structure formation and accelerated expansion. 
Furthermore, motivated by reasonable physical assumptions, we have been able to relate 
both dark components of the cosmological concordance model to the instability sectors 
of the dynamical system. 

As for the cases where the FLRW backgrounds are found to enjoy stability, we have 
to add a disclaimer. Despite their stable character, these FLRW backgrounds might not 
serve as reliable approximations in the role of physical backgrounds. As an example we 
look at the region II. b (green) in figure 1: we appreciate that the background experiences 
a transient acceleration period in between an initial and an asymptotic Friedmann-like 
state. Such a situation shows that FLRW backgrounds, even if globally stable, do 
a priori not provide a correct approximation for the physical background during the 
entire evolution. 

A last remark concerns the use of exact scaling solutions for the class of averaged 
models analyzed. Indeed, there is no proof that scaling solutions depict all the 
instabilities of FLRW backgrounds, but it is also highly improbable that this is the 
only class of averaged models where instability occurs. Since our point was to explicitly 
demonstrate the possible unstable character of Friedmann-like domains, the use of 
scaling solutions is perfectly legitimate, though not exhaustive. In other words, our 
study is restrictive in this sense, since scaling solutions, albeit exact, are only indicative 
for the existing instabilities. One could also speculate that scaling solutions are 
acceptable building blocks to describe any backreaction behavior, as any analytical 
function is the sum of a Laurent series; however the correspondence is not that trivial 
due to the nonlinearity of the problem. 
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